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Abstract 

We derive general expressions for the neutrino dispersion relation in a 
magnetized plasma with a wide range of temperatures, chemical poten- 
tials, and magnetic field strengths. If the electron and proton chemical 
potentials vanish, as in the early Universe, there is no magnetization 
contribution to the neutrino refractive index to leading order in the 
Fermi coupling constant, contrary to claims in the recent literature. 
Therefore, as long as the magnetic field satisfies B <T 2 , the neutrino 
refractive index in the early Universe is dominated by the standard 
"non-local term" . If neutrinos are Dirac particles with magnetic mo- 
ment jU, then their right-handed components are thermally populated 
before the nucleosynthesis epoch by magnetically induced spin oscilla- 
tions if fiBo > 10~ 6 /Ub gauss, where fiB = e/2m e is the Bohr magneton 
and -Bo is a large-scale primordial magnetic field at To ~ 1 MeV. For 
a typically expected random field distribution, even smaller values for 
IiBq would suffice to thermalize the right-handed Dirac components. 



I. INTRODUCTION 



If neutrinos carry magnetic or electric dipole or transition moments, they can 
spin-precess into other spin and/or flavour states in the presence of external magnetic 
fields. For example, if neutrinos were Dirac particles with a magnetic dipole moment 
/i, the active left-handed states could spin-precess into the otherwise sterile right- 
handed ones. It has been speculated that this effect can explain the deficiency of 
the measured solar neutrino fluxes, and it certainly can be important for supernova 
physics where large magnetic fields are known to exist 0. Further, it has been 
recognized for a long time that primordial magnetic fields of sufficient strength 
would couple right-handed Dirac neutrinos to the cosmic thermal heat bath and 
thus cause these "wrong- helicity" states to be thermally populated M. This effect 
would enhance the expansion rate of the Universe at the epoch of nucleosynthesis 
and thus modify the standard scenario of the formation of the light elements, in 
potential disagreement with the observationally inferred abundances. 

The original discussions of this cosmological effect || did not take into account 
neutrino dispersion, which at that time had received only marginal attention. Later 
on, it became clear that even though the neutrino dispersion relations in vacuum and 
in media are very close to that of massless particles, any deviation from the latter 
may cause significant modifications of spin or flavour-oscillation processes. A first 
assessment of medium-induced dispersion effects for early- Universe magnetic spin 
oscillations was provided in Ref. 0. In addition, however, one has to worry about 
neutrino collisions during the oscillation process. A formalism for the simultaneous 
treatment of oscillations and collisions was pioneered in Refs. @,[|, and was refined 
in terms of quantum-kinetic equations in Refs. 0. A quantum-kinetic treatment of 
the early-Universe magnetic oscillation problem was provided in a recent series of 
papers [TTHlQfl. 

Because even fine points of the neutrino dispersion relation are important for 
oscillation phenomena, one naturally wonders if the assumed presence of a strong 
magnetic field may cause a spin polarization of the electrons and positron in the 
medium, which in turn may act as a new contribution to the dispersion relation. 
Semikoz and Valle |8| claim that this is the case even for zero chemical potential, and 
that this effect dominates the neutrino dispersion relation for the physical conditions 
relevant in the early Universe. 

Upon closer inspection, however, we find that this dispersion relation is based 
on an unfortunate sign error. In a charge-symmetric plasma, the magnetization 
part of the local self-energy terms cancels between electrons and positrons rather 
than adding, as claimed by Semikoz and Valle ||. While the correct sign can be 
understood by a simple physical argument (Sect. |11B 2 ) and from the requirement 
of CPT invariance (Sect. [II B 3|) , we take this opportunity to provide the neutrino 
dispersion relation in a magnetized medium for arbitrary electron chemical potential 
and magnetic field strength. The correct sign is then a consequence of our completely 
general and formal derivation, which leaves no room for ambiguities. Our general 
expressions may also be of interest in the context of neutrino spin oscillations in 
supernovae, where strong fields and very degenerate electrons occur. Surprisingly, 
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we find that even for arbitrary field strengths our expressions are very similar to 
those derived by D'Olivo, Nieves, and Pal 



13 in the weak- field limit. 



Neutrino dispersion in a magnetized medium may be viewed from a somewhat 
different perspective where one considers an effective neutrino electromagnetic form 
factor, or vertex function, induced by the presence of the medium Jl3|^T6|. Various 
components of this vertex function, which is a Lorentz tensor, may be interpreted 
as certain effective neutrino electromagnetic multipole moments. In this language, 
neutrino dispersion in a magnetized medium is represented by a medium- induced 
effective neutrino magnetic dipole moment, which naturally leads to an energy shift 
in the presence of a magnetic field.Q The results of Refs. []T3"| , |i4]| imply that in 
a charge-symmetric plasma this dipole moment vanishes, in agreement with our 
present calculations and arguments. The same conclusion was reached in an early 
paper by Semikoz 
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in conflict with the later finding of Semikoz and Valle ||. 
In Sect. [TI| we derive general expressions for the neutrino dispersion relation in a 
magnetized medium, and we derive the relative sign of the magnetization effect by 
a direct physical argument. In Sect. [TIT] we investigate the efficiency of primordial 
neutrino spin oscillations in view of the correct neutrino dispersion relation in a 
magnetized plasma which, in the early Universe, is well approximated by the dis- 
persion relation of an unmagnetized medium. Section [IV] is devoted to a summary 
and discussion. 



II. NEUTRINO DISPERSION IN MAGNETIZED MEDIA 
A. General Self-Energy Diagrams 

In order to derive a general expression for the neutrino dispersion relation in a 
magnetized medium we observe that, to lowest order, the self-energy is given by 
the tadpole and bubble diagrams shown in Fig. [I]. To be specific we shall derive 
the dispersion relation for electron neutrinos; more general cases can be inferred by 
simple substitutions. 

The neutrino self-energy contribution from the tadpole diagram with an arbitrary 
fermion loop is 

-^tadpole = ~\ (^—] tl[i a {c V - C Alb )lS{x,x)} iD^^L , (1) 

4 VcostV/ 

where g is the weak gauge-coupling constant and 6w the weak mixing angle. We use 
the notation R = |(l+7s) and L = — 75). Further, D^A) is the Z-boson prop- 
agator while S(x, y) is the coordinate-space propagator for the background fermion. 



The use of an "effective magnetic dipole moment" to describe the neutrino energy shift 
in a magnetized medium is somewhat misleading, because the 7-structure of the vertex 
function is not that of a magnetic dipole interaction. Among other differences, only left- 
handed states experience any shift at all. 
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FIG. 1. The tadpole and the bubble diagrams. 



For a charged Dirac spin-| particle in the presence of an external magnetic field, 
S(x, x) is given in Appendix |A]. Our prime example is electrons for which the weak 
coupling constants are cy = — \ + 2 sin 2 9w and ca = — f ■ 

The bubble diagram contributes only for a background of charged leptons from 
the same family as the test neutrino. In our specific case of a test v e in the presence 
of an e + e~ plasma, we find 

- bubble = {^j Rj^i l a iS{v) iD%(k - p) j?L . (2) 

For a neutrino background from the same family as the test neutrino, there is a sim- 
ilar diagram with a Z—u-loop that can be obtained by replacing g 2 — > g 2 /2 cos 2 6^ 
and mw — ► ^z- 

The tadpole diagram provides only a local contribution, i.e. the gauge-boson 
propagator is taken at the energy-momentum transfer A = so that we could have 
used an effective low-energy four-fermion interaction. The bubble diagram, however, 
involves the gauge-boson propagator at a non-vanishing A so that there is a non- 
local term in the self-energy. Even in extreme astrophysical sites, such as neutron 
stars, the relevant energies are so low, and the chemical potential so high relative to 
the temperature, that the bubble diagram is dominated by the local term. However, 
the local term vanishes identically in a charge-symmetric plasma. Therefore, in the 
early Universe the neutrino self-energy is dominated by the non-local part of the 



bubble diagram [11 . 



B. Local Terms 
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1. Formal Derivation 



In order to derive the electron-neutrino dispersion relation in a magnetized 
medium explicitly, we begin with the local contributions. To this end we expand 
the gauge-boson propagators in powers of the energy-momentum transfer A, 

gg z (A) = + ** A ' 7 A ° A " + O (*) . (3) 

m w,z m w,z \ m w,z ) 

The first term, which is the only one contributing to the tadpole, gives the local 
part of the self-energy. 

Using the charged-fermion propagator in an external magnetic field (see Ap- 
pendix |A| for more details), the tadpole yields for a plasma consisting of electrons, 
neutrinos, and nucleons: 



-N n . n + 2 £ N^._ p . - (1 - 4 sin 2 6 W ) (iV e _ e - - N p . p ) 



7o 



+ N e _ s B-j \L , (4) 



_Gf( 

J tadpole v/2 I 



where S is a unit vector in the external 5-field direction. Further, iVy_j denotes the 
net number density of fermions /, i.e. the total number density of fermions / minus 
that of antifermions /. For neutrinos, only the number density of left-handed states 
(superscript L) is counted, which is identical to the total number density unless the 
right-handed degrees of freedom have been populated by, say, magnetically induced 
spin oscillations. Usually, the standard electron and proton terms cancel against 
each other in a charge-neutral plasma where N e _ E — N p _p = 0. 

In the magnetic tadpole term, N®_ s is the net number density of electrons in the 
lowest Landau level. Of course, the exact cancellation of all higher Landau levels ap- 
plies only to Dirac fermions which do not carry anomalous magnetic dipole moments. 
This approximation is not justified for nucleons, which carry large anomalous mag- 
netic moments so that their polarization does not cancel between the higher Landau 
levels which are not degenerate. However, unless the field is extremely strong or 
the temperature much higher than the nucleon masses, the nucleon magnetization 
is suppressed by their heavier masses relative to electrons. Because in the present 
paper we are primarily interested in early-Universe physics between the QCD phase 
transition and Big-Bang nucleosynthesis (BBN), nucleons can certainly be ignored 
with regard to neutrino dispersion effects. 

In addition we need to consider the bubble diagram, which yields a local contri- 
bution from electrons and electron neutrinos of 
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iV* + N e _ s ) 7o - N° e _ s B 7 L . (5) 



^bubble — ~\/2 ^ 

Nucleons never contribute to this term. 

In the weak-field limit our results agree with those found in Ref . [ D^] , except for 



the overall sign which is related to the convention in Ref. [13] that e < for electrons. 
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The approach in Ref. |13j was strictly perturbative in that a plane-wave basis for 
the fermions was used instead of Landau levels. We stress that exact expressions for 
quantities such as the magnetization or the magnetic susceptibility do not in general 
admit a power-series expansion in B, forcing one to use Landau levels as external 



states However, when a quantity does admit a power series expansion, it is not 
too surprising that the linear term of the exact result agrees with a perturbative 
calculation based on plane-wave states. 

Our magnetic neutrino self-energy terms apply for B <C m^, but B may well be 
large compared with other scales in the problem, such as the electron mass or the 
temperature. Even for such large fields the linear term actually gives the complete 
result. This surprising finding is traced to the fact that only the lowest Landau level 
contributes and that N®_ s is strictly linear in B. It must be noted, however, that the 
presence of the field affects the phase-space distribution of the charged fermions and 
thus the relationship between chemical potential and density. Therefore, one must 
specify if the charged-particle densities or their chemical potentials are held fixed in 
order to specify the functional dependence of the neutrino dispersion relation on B. 

The dispersion relation for left-handed electron neutrinos in a magnetized plasma 
is obtained by taking the determinant of jk — S tadpo i e — S bubble . We find 

E± = ±k = ±a+ \k-b\ , (6) 

where ± refers to v e and i> e , respectively. Further, 

I/vr _ + V N L + N L - + N - 

2 ly n—n i / , J-^ Ui - Ui ~ -*V e — i> e ' Jv e— e 
i=e,/j,,T 

-{\-2sm 2 6 w ){N e _- e -N p _ p ) , 
6 " \Nl- e B . (7) 



V2G F 



V2G 



It is the medium- and field-induced breaking of Lorentz invariance that generates 
a non-trivial dispersion relation, or refractive index, for neutrino propagation. In 
a charge-neutral plasma the term proportional to (N e _ s — N p _p) vanishes. Again, 
there is a small nucleon contribution to b which we have neglected. We stress that 
it is a slight abuse of language to call b magnetization because only the spin part of 
the magnetization enters, not the orbital part. Note further that the spin is not a 
conserved quantity and only the lowest Landau level is a spin eigenstate. 



2. Physical Derivation 

Because the local magnetization contribution to the refractive index is contro- 
versial in the literature, it is useful to provide a more physical derivation where the 
absolute sign, and the relative sign between the electron and positron terms, become 
more directly apparent. To this end we may start directly from the four-fermion 
neutrino vertex with a charged lepton £: 

n int = V2G F ^ ula L^ u ^a a (gv-gAi 5 )^i • (8) 
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Here, the effective weak neutral-current coupling constants gv, A are identical with 
cy y A unless £ is from the same family as the neutrino, in which case gy t A — cv,a + 1 
because the Fierz-transformed charged-current mimics a neutral-current interaction. 

The neutrino self-energy is found by calculating the expectation value (\I/£7 a (gy— 
gA^^e) m a background bath of fermions £. In an unpolarized, isotropic medium 
only the zeroth component of the vector current contributes and yields the standard 
result. A magnetically induced polarization of the charged background fermions, 
however, causes the axial current to obtain a non-vanishing expectation value. 

For ultrarelativistic charged fermions the expectation value of the chirality op- 
erator 75 is identical with that of sign(g)p ■ BX, an observation that establishes a 
simple relation between chirality and the magnetic quantum number A of the Lan- 
dau levels. It implies that the axial-vector contributions cancel between charged 
fermions with the same momentum but opposite A. The Landau-level energies 
E n,\,p* = m 2 + p 2 z + \qB\(2n + 1 - A), with n = 0, 1, 2, ... and A = ±1, are de- 
generate between the levels (n, A = +1) and (n — 1, A = —1) except for the lowest 
level (n — 0, A = +1), which is not matched by a lower level with opposite mag- 
netic quantum number. Therefore, only the lowest Landau level contributes to the 
expectation value of the axial- vector current. 

A negatively charged ultrarelativistic £ in the lowest Landau level, moving along 
the -B-field, has its magnetic moment parallel to B, a spin opposite to B, and 
therefore negative chirality. Hence for such a state 

i^a{gv - gAjsW = (^n^e)(gv + 9a) = B( gv + g A ) . (9) 

If £ moves in the opposite direction we get —B(g v — gA) so that the vector part 
averages to zero for each momentum mode separately if the phase-space distribution 
is reflection-symmetric along B. An anti-£ (£) moving along the .B-field has its 
magnetic moment also pointing parallel to B, but its spin in the opposite direction 
relative to an I with the same momentum along the field, since the charge is opposite. 
Therefore, the helicity of £ is opposite to that of £ and thus their chiralities are equal. 
The expectation value corresponding to Eq. is then —B(g v + gA)- Similarly we 
get B(gy — gA) for an £ moving in the opposite direction. Multiplying with the net 
number density of €s and £'s in the lowest Landau level we obtain 

X v = -g A y/2G F N°_ l B->yL , (10) 

where the final minus sign comes from the contraction of space-like indices in Eq. (§) . 
While our simple derivation was based on the notion of ultrarelativistic charged 
leptons, this result holds true even for non-relativistic ones as follows from the 
formal derivation in the previous section. 



The absolute sign of the energy shift, which differs from the one found in Ref. fL3 
can be checked by comparing Eqs. (f|) and (|J) with Eq. fllPl) . If the leptons are of a 
family opther than the neutrinos, gy t A — cv,a, and one must compare Eq. (|T0| ) with 
the tadpole term alone. For leptons of the same family, g V A = cy,A + 1, and the 
sum of the bubble and the tadpole diagram should be compared to Eq. ([TTJ]) . The 
relative sign between the electron and positron terms also follows directly from this 
simple derivation without ambiguity. 
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In the early Universe, the numbers of particles and antiparticles are believed to 
be identical to within about 10 -9 . Therefore, the plasma was effectively charge- 
symmetric. Our general results, Eqs. (|]) and (||), reveal that, to leading order in 
m^ 2 , there is no magnetization contribution to the neutrino refractive index in 
such an environment, contrary to what has been claimed by Semikoz and Valle || 
who found that the fermion and antifermion terms in Eqs. @ and (||) add rather 
than subtract. It is correct as in Eq. (2.9) of Ref. M to identify the relevant spa- 
tial part of the axial current (^/ e / j , ~f^ e ) with the difference between the electron 
and positron magnetizations, but in the manipulations leading from their Eq. (3.3) 
to (3.6) Semikoz and Valle have unfortunately picked up an incorrect minus sign. 
Therefore, at epochs before nucleosynthesis the non-local neutrino refractive terms 
remain more significant than the local ones JTTJ, even in the presence of strong 
magnetic fields. 



3. CPT Argument 

The vanishing of the local contribution to the neutrino self-energy in a CP- 
symmetric plasma can also be deduced from a direct symmetry argument. To this 
end we assume that the background plasma is CP symmetric, and in addition we 
assume that it is in a stationary state so that it is also symmetric under the time- 
reversal operation T. Since CPT is strictly conserved in our theory, and the magnetic 
field is CPT invariant, it follows that neutrinos and antineutrinos of a given momen- 
tum must experience the same medium-induced energy shift, i.e. their self-energy in 
the medium must be the same. Put another way, the expectation value of ty v Y& v 
must be the same for neutrino and antineutrino states of equal momenta. 

At one-loop level, the general form of the self-energy operator £ in a magnetized 
medium is [[13] 



£ = R(aW + bu^ + cB^L . (11) 

Here, k is the four-momentum of the test (anti)neutrino, u is the four-velocity of 
the background medium, and = ^e^ va/ 3U u F al3 is a covariant expression for the 
external electromagnetic field which is a pure B-field in the rest frame of the medium. 
The coefficients a, b, and c are functions of the scalars k 2 , uo = k ■ u, and k ■ B. 

Under CPT the current ^ u , y fl ^ l , and the four-momentum k change sign. (Recall 
that the Dirac eq. implies ^ ' Jz^^i^ ' v = m^ u ty u , and that is invariant under 

CPT.) However, the four- vectors u and B are invariant under CPT. It is important 
to observe here that u is not an operator for ^> v since it is just fixing the new 
reference frame. From Eqs. ([I]) and (0) it is clear that the local contribution to E 
is independent of k so that the coefficient a must be zero while b and c must be 
constants. However, because and B^ are even under CPT while {^ u ^^! v ) is 
odd, and because (^/ U T,^/ U ) is required to be even, we find that b and c must be 
zero. The coefficient c is related to a medium-induced effective neutrino magnetic 



dipole moment. In Ref. [20] it was already shown on the basis of the same argument 



that such a dipole moment must vanish. Thanks to this argument no contributions 
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to Siocai can arise from strong-field corrections to the W propagator in a CPT 
symmetric plasma. Non-local terms which are odd functions of scalars that are odd 
under CPT, namely u and k ■ B, are not required to vanish. 



C. Non-local Terms 



In a charge-symmetric plasma, all of the local self-energy terms given in Eq. (|7|) 
vanish so that the second term in the expansion of the gauge-boson propagator in 
Eq. (H) dominates. We shall concentrate on the case where m e < T < m w and 
B < T 2 . In the early Universe, these are quite reasonable approximations between 
the QCD phase transition and nucleosynthesis. Repeating the calculations for the 
bubble diagram we then obtain 



J bubble 



(k) 



7V2vr 2 G F T 4 



45 m 2 z 



1 + 



2m 2 z 



loh - ~.lk L 
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V2G F T 2 
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w 



7ofco + (B • j)(B ■ k) 



(12) 



where cr is a vector of Dirac spin matrices defined by B ■ cr = jF^^^j^], where 
F^ u is the field strength tensor. The resulting dispersion relation takes the form 



E 4 



±kn 



7V2tt 2 G f T 4 
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(13) 



The first part agrees with the result of Ref. [llj; it is the same for v e and v e . The 
i?-dependent energy shift is anisotropic and opposite for v e and v e . However, it 
remains subdominant compared to the isotropic term as long as B < T 2 . 



III. SPIN OSCILLATION IN THE EARLY UNIVERSE 

A. Neutrino Depolarization Rate 

As an application of our results we consider explicitly the case of Dirac neutrinos 
with a magnetic moment \i. In the presence of a primordial magnetic field the 
thermally populated left-handed (l.h.) states can spin-precess into the otherwise 
sterile right-handed (r.h.) ones, thus which will be populated as well. This process 
of populating the "wrong-helicity" neutrino states is treated theoretically by virtue of 
a Boltzmann-type kinetic equation, which includes neutrino oscillations as discussed 
in Refs. 0. However, for a simple estimate one may ignore the detailed evolution 
of the individual momentum modes and rather study an average evolution of the 
overall spin-polarization vector P of the entire ensemble. In this simplified approach 
the global spin-polarization vector evolves as |J 



S 



d t P = V x P - DP T , 



(14) 



where V is a vector of effective magnetic interaction energies, D a damping rate 
due to collisions, and Pt the "transverse" part of the spin-polarization vector to be 
discussed below. 

In the absence of a medium, the damping rate vanishes and the effective interac- 
tion energy for ultrarelativistic neutrinos is V = 2p,B^, where Bt is the component 
of the B field which is transverse to the neutrino direction of motion ■ Be- 
cause only the transverse magnetic field matters, in vacuum the neutrino helicity 
can be reversed entirely by spin precessions. Put another way, l.h. and r.h. states 
are maximally mixed by the presence of a magnetic field, independently of the field 
direction with respect to the neutrino direction of motion, unless B^ vanishes ex- 
actly. Of course, the precession time depends on the magnitude of B^ and thus on 
the relative field direction. 

The first impact of a medium is that it endows the active (l.h.) neutrino states 
with a nontrivial dispersion relation, while the sterile (r.h.) ones remain unaffected. 
Because the particle-antiparticle asymmetry in the early Universe is thought to be 
of order 10~ 9 for all species, i.e. small, the dominant contribution to the neutrino 
refractive index is the non-local term that was first identified in Ref. [11]. The 



discussion in Sect. [11 H reveals that even in a magnetized charge-symmetric plasma, 



the additional neutrino refractive term from the B field is rather small so that 
the standard isotropic term continues to dominate, in agreement with the treatment 
of Ref. @. We expect this to remain true, and that the first term of Eq. fll2|) remains 
good as an approximation even for temperatures not much higher than the electron 
mass which are relevant at the time of the BBN. However, in Ref. the impact 
of the damping term was not properly discussed. The later systematic studies of 
kinetic equations for oscillating neutrinos were not available at that time. 

In order to identify V and D relevant for the conditions of the early Universe we 
begin with the energy difference between l.h. and r.h. neutrinos of flavour i — e, fi, 
or r which is i?i.h. — £7 r .h. = — Here, E is the unperturbed energy, which agrees 
with E T _Yi. because r.h. neutrinos do not experience any energy shift in the medium. 
Assuming the mass of the neutrino to be much smaller than the temperature, it is 
easy to extract the coefficient £ from Eq. flI3|): 

_ 8^/2 G F ( Pv l+Dl p e+ 



Tflr 




+ ■ (15) 



Here p^.+p, is the energy density in l.h. neutrinos plus antineutrinos of flavour £ 
while p i+ i is the energy density in the ^-flavoured charged leptons plus antileptons. 
The coefficient £ has the same sign for neutrinos and antineutrinos as test particles. 
For t neutrinos in the early Universe it is dominated by the P„ t+9t term because the 
presence of r leptons is suppressed by a Boltzmann factor e~ mr//T . 

In a magnetized charge-symmetric plasma the spin-polarization vector of an 
ultrarelativistic neutrino or antineutrino of energy E evolves according to Eq. (|T4"D 
with 
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\V L \=ZE , (16) 

where T and L are understood to be transverse and longitudinal relative to the neu- 
trino direction of motion. Since we are using neutrino helicity states, the direction 
of spin-quantization is identical with the direction of motion. Therefore, the tilt of 
V relative to the direction of motion is twice the effective in-medium mixing angle 
H between l.h. and r.h. states: 

tan2fl = ^ = ^ , (17) 

where Vt,l — |^t,l| and = \Bt\- Thus, for sufficiently weak magnetic fields 
the l.h. and r.h. states are effectively de-mixed so that l.h. states spin-precess only 
partially into r.h. ones. Put another way, the spin precession is about the direction of 
an effective magnetic field B c q = V/fi which is no longer transverse to the direction 
of motion. 

The second effect of a medium is that l.h. neutrinos scatter, thereby interrupting 
the precession process. A collision essentially amounts to a "measurement" of the 
helicity content of a given neutrino because the l.h. component is scattered out of 
its previous direction of propagation while the r.h. component moves on unscathed. 
This implies that every collision resets the neutrino into a helicity eigenstate and the 
oscillation process begins from scratch. Collisions thus destroy the phase coherence 
between the l.h. and r.h. component of a neutrino state, which amounts to a damping 
of the transverse part Pt of the polarization vector. In the present situation where 
the r.h. component does not interact at all, the damping rate D is found to be half 
the collision rate of the l.h. component |||6| so that, in the early Universe: 

D = h^GlT*E . (18) 

Here, /d is a numerical factor, which was found to be unity for fi- or -/"-flavoured 
(anti) neutrinos in a background medium of e and all sequential (anti) neutrinos 
]IIJ. Corrections from the magnetic field are expected to be small for eB < T 2 . 

For an estimate of the rate of depolarization Tdepoi of the initially l.h. (anti)neu- 
trino population, we turn to thermal averages of the refractive and damping terms. 
The average energy is (E) « 3T for a given neutrino species where the equality would 
be exact if the neutrinos would follow a Maxwell-Boltzmann distribution instead of 
a Fermi-Dirac one. Then 

(D) = h 7 -^G 2 F T 5 (19) 

is the average damping term. For the refraction term in Eq. (|15|) we note that the 
energy density in one flavour vi of l.h. (anti)neutrinos is p Ue +u e = (77r 2 /120) T 4 . 
Further, m^ 2 = (\f2G-p /na) sin 2 W cos 2 6 W where a ~ 1/137 is the fine-structure 
constant. We will always approximate sin 2 6w = 1/4. With (E) = 3T we thus find 

(Vi)^/ L ^G 2 T 5 , (20) 
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where /l = 1 + (Pi+i/ Pu e +p e ) ( m z/mw) 2 is a factor to account for the possible 
presence of charged leptons of flavour £, which would also contribute to Eq. (pf). 
For p and r neutrinos the lepton densities are small and we have /l = 1, while for 
e neutrinos /l ~ 3.6. 

These results are enough to determine that the evolution Eq. fll4| ) of the neutrino 
polarization vector is weakly damped, i.e. that it typically precesses several times 
between collisions. The oscillation or spin-precession frequency is identical to V = 
(V T + Vl) 1 ' 2 > V L . With Eqs. (HD and © we find 

O) h 5« (21) 



(Vl) h 2 

The average period of spin precession is approximately 2tt/ (Vl,) so that there are at 
least about 10 revolutions between collisions. 

In order to understand the solution of Eq. (|T4] ) in the weak-damping limit we 
multiply both sides with P/P 2 , which leads to d t P/P = -D (Pt/P) 2 with P = \P\ 
and Pt = \Pt\- I n the weak-damping limit we may use a precession-averaged Pt, 
which is found by taking the transverse part of the projection of P on V. Elementary 
geometry yields Pt/P = cos 29 sin 29 so that in the weak-damping limit 

d t P/P = -D cos 2 26 sin 2 29 , (22) 



where the effective mixing angle is given by Eq. (|TT|). Assuming that it is small, so 
that cos 29 pa 1 and sin 29 pa tan 29 we thus find 



(2fiB T ) 2 (D) f D 400 a 2 fi 2 B 2 



depol ~ /T/2 \ ~ f 2 



(V 2 ) f 2 7n GpT 5 



(23) 



for the average depolarization rate r depol = (D cos 2 29 sin 2 29). We have freely factor- 
ized the thermal averaging process, and we have assumed a homogeneous magnetic 
field. 



B. Comparison with Expansion Rate 

In order to decide whether the depolarization of the initially l.h. neutrino ensem- 
ble is ever complete during the cosmic evolution, we need to compare Tdepoi with 
the expansion rate H. If at some epoch r depol > H then r.h. neutrinos have approx- 
imately reached thermal equilibrium at that time. If this epoch falls between the 
QCD phase transition at Tq CD ph 150 MeV and nucleosynthesis at T B bn ~ 1 MeV, 
then a significant impact on the primordial light-element abundances would have to 
be expected. 

According to the Friedmann equation the expansion rate is given by H 2 = 
(87r/3) p/rripY, where p is the energy density of the Universe at a given epoch and 
m-pi = 1.22 x 10 19 GeV is the Planck mass. In the radiation-dominated epoch, the 
energy density is p = g (7r 2 /30) T 4 , with g the effective number of thermally excited 
degrees of freedom. Between the QCD and BBN epochs we need to count photons, 
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e^, and the l.h. sequential (anti)neutrinos so that g = 43/4. Therefore, we need to 
require that at some epoch r depol exceeds 

J ff = ^(437T 3 /45) 1 / 2 r 2 /m P1 , (24) 

where f H = (^4/43) 1 / 2 . The energy density of the magnetic field should be added 
to p, but since it is at most of the same order of magnitude as the radiation energy 
it can be absorbed into /# without changing our analysis significantly 

In order to perform this comparison we need to understand the scaling law of an 
assumed magnetic field distribution under the cosmic expansion. Flux conservation 
indicates that B oc Rr 2 or B oc T 2 . This would be the complete scaling if the 
magnetic field were homogeneous. In practice, any primordial field distribution is 
expected to be complicated and noisy, so that the effective B\ in Eq. must 
be interpreted as a suitable average. Essentially, each neutrino oscillation process 
"measures" the magnetic field linearly averaged over distances corresponding to 
the oscillation length so that we need to estimate (-Bt)l osc = \ Il osc \d x \B^(x), 
where the integral is over a neutrino oscillation path. This linear average is relevant 
because the spin-precession equation (Eq. (0)) is linear in B^. Further, one must 
take an ensemble average over all oscillation paths at a given epoch. Here, the 
average should be taken over the quantity (-E?t)| osc since the expression for the 
depolarization rate in Eq. ( p3|) is quadratic in B^. The co-moving oscillation length 
L osc increases with time so that at later times the effective field strength is averaged 
over larger co-moving domains, reducing the effective field strength below the naive 
T 2 scaling law. 

As a simple model for the 5-field scaling we assume a power law in co-moving 
coordinates of the form []nj 

«B r »^«B T >i> 1 / 2 = J B i (|) a ^ 7 , (25) 

where T is the temperature at a reference epoch which we take to be the BBN time 
with T Q — 1 MeV, d is the co-moving size of a typical domain over which the field 
is correlated, and is the field strength in such a domain. The average ((. . 
indicates a linear averaging over the length scale L and a root-mean-square average 
over all oscillation paths. The exact I?-field scaling depends on the mechanism 
of initial generation and the evolution of the complicated magnetohydrodynamic 



equations ||26|| . Therefore, the power law in Eq. (|25| ) should only be taken as a toy 
model for L ^> d. In particular, the merging of two domains would change both d 
and Bd while we take the combination B^T^d 1 to be constant here. 

Because oscillating neutrinos measure the magnetic field over an oscillation 
length scale, it is natural to use L = L osc , which is of order H^ 1 at BBN . There- 
fore, we define a horizon-scale magnetic field at BBN by B = ((B^)) Losc taken at 
T = Tq. The condition Tdepoi ~ H for which r.h. neutrinos are certain to reach 
thermal equilibrium at some epoch T then translates into 

n ^ffnflY /2 7T (21077r\ 1 / 4 G F T 3 / 2 T 2 f L osc (T)V , s 
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The temperature dependence of the co-moving oscillation length can be determined 
from IV]- 1 < V^ 1 in Eq. © so that L OSC (T)/L OSC (T ) = (T /T) 4 . If we focus on 
the period between the QCD phase transition and BBN we have fn = 1, and the 
numerical coefficient in Eq. ( p6l) is 0.55 ~ 1/2. Then 

VB >^^T^T^ , (27) 
2a mpj 

where in addition we have used /d = /l = 1 appropriate for and v T . 

Evidently we need to distinguish between two generic cases depending on whether 
7 < 3/8 or 7 > 3/8. Beginning with the former, which would be applicable for a 
homogeneous field, the condition in Eq. (^7|) should be imposed at as low a temper- 
ature as possible in order to find the smallest necessary fiB sufficient to populate 
the r.h. states. Our entire discussion makes sense only as long as neutrinos scatter 
efficiently by weak interactions. At later times they may still spin-precess in the 
cosmic magnetic field, but this would have no further impact on the expansion rate 
of the Universe as the only effect would be to redistribute the total neutrino energy 
density between r.h. and l.h. states. Neutrinos freeze out at about T = 1 MeV, 
just before the BBN epoch. With T = T = 1 MeV r.h. neutrinos reach thermal 
equilibrium before BBN if 

Q rpJ /2 

I^B > F «7x 1(T 15 eV « 1.2 x 10~Vb gauss . (28) 
2a m P1 

This requirement is essentially identical to what was found in Ref. 0, even 
though the interplay between oscillations and collisions was not treated there. It was 
demanded that the mixing angle should be large, and that the damping rate should 
be small compared with the spin-precession rate, conditions which are sufficient, 
but not necessary to achieve thermal equilibrium. Here we found that we are always 
in the weak-damping case. If we take damping effects into account according to 
Eq. (|T4]), the required magnitude for fiB at the critical epoch around neutrino 
freeze-out implies that the mixing angle is not small. Therefore, either treatment 
leads to roughly the same answer. The underlying reason for this coincidence is 
that in the early Universe the dispersive and the absorptive parts of the neutrino 
refractive index are of the same general magnitude, i.e. they are both second order 
in Gp. Then, at the critical epoch around neutrino freeze-out, the time scales (VL), 
(D), and H are all about the same to within numerical factors. 

The assumption that the magnetic field is a slowly varying function (7 « 0) on 
the scale of the Hubble radius at nucleosynthesis is not physically very likely. In fact, 
an ubiquitous mean field would be incompatible with the observed cosmic isotropy 



if its present strength were larger than about 10 gauss |24]]. Furthermore, it is 
a general feature of models predicting magnetic field generation during primordial 
phase transitions [25] to forecast random magnetic fields at the end of the transition, 



in domains having a typical size d <C H~ l . Although magnetohydrodynamical f26 
and dissipative effects |27]] can cause the ratio d/H" 1 to grow during the cosmic 
expansion it may still be much smaller than unity at the BBN time. At that epoch, 
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the neutrino oscillation length is not much smaller than H~ l , so that the neutrino 
probes a number of field inversions before one spin precession is complete. 

If the magnetic field performs a random walk along each neutrino trajectory, the 
average transverse field decreases with the square root of the length scale. Therefore, 
one would expect that 7 = 1/2, whence it appears more natural that 7 > 3/8. 

In this second generic case the condition Eq. (p7|) is easiest to fulfil at early 
times. Typically, the earliest useful epoch is just after the QCD phase transition 
at T ~ 150 MeV. Then, because 7 > 3/8 by assumption, the required value for 
fiB will be smaller than Eq. ( [2"8"D by a factor (T /T) 47_3//2 , which for 7 = 1/2 is an 
order of magnitude. Therefore, Eq. (|28| ) is a conservative requirement in the sense 
that this value for [iBq is certainly sufficient to populate r.h. neutrinos before BBN, 
but a smaller value may suffice, depending on the exact scaling law of the effective 
5-field. 

In our derivation we have assumed that the effective mixing angle is small, a 
condition that we now need to verify. From Eq. ( [P7D we need to demand that 
2(fiB T )/(V L ) < 1 or 

»B Q < h ^~ GlT^T 2 ^ . (29) 
80 a 

This condition is most difficult to fulfil at late times, unless 7 > 3/4. Therefore, 
it is enough to check it at T — Tq — 1 MeV. At that temperature it amounts 
to [iB < 5 x 10~ 15 eV. A comparison with Eqs. (|28| ) and (|29|) reveals that our 



assumption of a small mixing angle has been marginally consistent for 7 < 3/8, 
and safe for 3/8 < 7 < 3/4. Assuming that the mixing angle is large amounts to 
ignoring refractive effects. This leads to a requirement similar to Eq. (|28| ) for r.h. 
neutrinos to reach thermal equilibrium. 

The magnetically induced spin-oscillation of neutrinos in the early Universe has 
been discussed in several recent papers 0. While some of them discuss the im- 
portance of neutrino refractive effects at length, this effect does not always matter 
in their final result. The difference to our treatment is that we study the effect of 
correlated domains with finite sizes while these papers use the limit where the fields 
in different points are uncorrelated, (Bi(x)Bj(y)) ~ 5ij5^(x — y). In that limit, the 
magnetic field is assumed to consist of very small domains of random magnetic field 
strength and direction. Therefore, the main difference between our discussion and 
that of Refs. consists of the assumptions about the magnetic field distribution, 
and the kinetic treatment adequate for those assumptions. 



IV. DISCUSSION AND SUMMARY 

We have studied magnetically induced spin precessions of Dirac neutrinos in the 
early Universe. To this end we have derived expressions for the neutrino dispersion 
relations in magnetized media which are valid for field strengths B up to about rn^r. 
In the weak-field limit, our results agree with those of D'Olivo, Nieves, and Pal ITS" 



apart from an overall sign. In a charge-symmetric plasma, there is no magnetization 



contribution to the neutrino refractive index to lowest order in mJ?, contrary to the 
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claim of Semikoz and Valle [Q. Besides a formal derivation, we have shown how to 
obtain the magnetic refraction term in a direct and simple physical fashion, which 
establishes without ambiguity the absolute sign, and the relative sign between the 
electron and positron contributions. 

Our analysis indicates that r.h. Dirac neutrinos would be thermally populated by 
spin oscillations if [aBq > 10~ 6 /zb gauss, where \i is the assumed neutrino magnetic 
dipole moment, fis = e/2m e is the Bohr magneton, and Bq a horizon-scale magnetic 
field at T = 1 MeV, i.e. just before the epoch of nucleosynthesis. Depending on the 
spatial magnetic-field distribution on smaller scales, i.e. with sufficient power in 
smaller-scale field modes, even a smaller value of fiB would suffice to thermalize 
the r.h. states. 

In principle, r.h. neutrinos could also be populated by direct spin-flip collisions 
on charged particles or from annihilation processes involving virtual photons |p8 



The dipole moment needed to achieve thermal equilibrium for the r.h. states is 
fi > 0.5 x 10 -10 /!b- If the neutrino mass is smaller than 1 MeV, as we assume in the 
present paper, stellar-evolution bounds on neutrino dipole or transition moments 
are fi ^ 3 x 10 _12 /iB |Pl,p9|, so that the scattering mechanism cannot be effective in 
the early Universe. 

In the particle-physics standard model, neutrinos have no magnetic dipole mo- 
ments. However, if neutrinos have a Dirac mass m they automatically have a dipole 
moment fi/fiB = 3.2 x 10~ 19 m/eV. In other extensions of the standard model much 
larger values can be obtained. If one of the neutrinos would saturate the stellar- 
evolution limit, a primordial field at nucleosynthesis B ps 3 x 10 5 gauss would be 
enough to populate the r.h. degrees of freedom, and an even smaller field could 
suffice, depending on its spatial distribution. 

Unfortunately, direct observations of primordial magnetic fields are still lacking, 
although it was recently suggested that they may be detectable by observing their 
inprint on the cosmic rays or on the cosmic microwave radiation |31| . However, 



it may be useful to consider some recent hypotheses about the genesis and evolution 
of primordial magnetic fields. Many of these propositions are motivated by the 
desire to explain the observed galactic and intergalactic magnetic fields as relics 
of a primordial cosmological field. Field strengths of order 10 -6 Gauss are a quite 
general character of the interstellar medium. Remarkably, this strength corresponds 
to an energy density equal to that of the cosmic microwave background radiation. 
Kronberg [32] suggests that this feature may be the result of an early equipartition 
between magnetic fields and radiation, a hypothesis that may have found some 



theoretical support (e.g. Ref. [j3~3~| ). If this were the case we could expect B « 
10 13 Gauss, a value which is not in contradiction with primordial nucleosynthesis 



considerations | 53}| . If such large fields were produced before nucleosynthesis, our 
result implies that even a dipole moment as small as about 10 _19 /xb would be enough 
to thermalize r.h. neutrinos. Thus, neutrinos with cosmologically interesting Dirac 
masses in the eV range would have sufficiently large dipole moments without further 
extensions of the standard model. 

It has frequently been argued that the observationally inferred primordial light- 
element abundances exclude significant novel contributions to the cosmic expansion 
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rate of the Universe at the nucleosynthesis epoch. At the present time, however, 
new questions concerning the reliability of the previously inferred abundances of 
deuterium and 3 He have arisen, and the overall consistency of BBN with all of the 
observations is not assured ||35|| . Therefore, at the present time one cannot assume 
that the observationally inferred primordial light-element abundances truly exclude 
one additional thermally excited neutrino degree of freedom at the nucleosynthesis 
epoch. Therefore, it is not the ambition of our present study to claim a new exclusion 
range for fiB , but rather to illuminate some of the important physical ingredients 
needed to understand magnetically induced neutrino spin oscillations in the early 
Universe. 
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NOTE ADDED 



Before circulating the present paper as an E-print we made a draft version avail- 
able to Drs. J.W.F. Valle and V. Semikoz who subsequently agreed that our expres- 
sion for the magnetically induced refractive index was the correct one. As a formal 
response they have now circulated a corrected version of their derivation [36| which 
explicitly confirms our finding. 



APPENDIX A: CHARGED-FERMION PROPAGATOR 



In order to calculate the tadpole and bubble diagrams in Sect. [0], we need an 
explicit expression for the electron propagator S(x', x") in the presence of an exter- 
nal magnetic field and an electronic plasma at non-zero temperature and density. 
We shall use two different methods here: Furry's picture for the local term and 
Schwinger's proper-time method for the non-local one. They give the same result 
for the local terms, but the Furry-picture result is more direct to interpret phys- 
ically. For the non-local terms it would be considerably more difficult to use the 
Furry picture. 

By the Furry picture we mean that the propagator is constructed explicitly as a 
sum over solutions to the Dirac equation in a given gauge. For a fermion with mass 
m and charge q (the electron having a negative charge q = — e < 0) the propagator 
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has been constructed in Refs. p2[^3j. For a magnetic field in the positive z-direction, 
in the gauge = (0, 0, — Bx, 0), it is given by 



iS(x, x) 



f+°° dp 




dp z 


i 


Loo ~2n J 


— oo 2n J 


-oo 2n 


U - Et Pz 



- 2«5{pl - ElJf F (p Q ) 



x I (Polo - Pzlz + m) a + I l}l (x,p y ) + a^I^i^i(x,p y ) 



2l\qB\ 'j + Ii ) i- 1 (x,p y ) -J-Il-.ljfaPy) 



(Al) 



whereQ 7± = ^[jx ± S ig n (<7-B)nvl an d er± = |[1 ± siga(qB)a z ]. Note that a z is 
understood to mean the Dirac spin matrix f [7x;7j/]- The Landau levels are labelled 
by I and their energies are Ef = m 2 +p 2 z + 2\qB\l. Further, /f(Po) = fpipo) @(Po) + 
/f (Po) Q(-Po), where /p(po) = (e 



3 ±(p -M)/T . 



are the usual occupation numbers 
for particles and antiparticles of a Fermi-Dirac distribution at temperature T and 
chemical potential \i. We have also used Ik,i(x,p y ) = Ik(x,p y )Ii(x,p y ) with 



h{x,p y ) 



■\qB\ 



1/4 



7T 



exp 



\qB\ 



x — 



Py_ 

qB / 



H, 



(A2) 



where Hi is a Hermite polynomial. In the lowest Landau level we define I_i = for 
consistency 

The dp y integration can be performed by using the completeness relation 



dpyI k (X,Py)ll(X,Py) = \qB\8, 



kl 



(A3) 



which also removes the x-dependence from the r.h.s. of Eq. ( |A1| ). In the end we are 
only interested in the thermal part, coming from the 5- function in Eq. (|Ai"|), so we 
drop the vacuum contribution from now on. After the dp z integration has been done 
using the 5-function we find 



iS(x, x) 



\qB\ 
4tt 2 



dp fF(Po) (loPo + m 

,2\ 



X 



Qipl-m 2 ) ~ Q(p 2 - m 2 - 2\qB\l) 

Pq — mr i=i 



Pq — m 2 



2\qB\l 



(A4) 



The appearance of the projection operator a + in Eq. (|A4j) is related to the fact that 
there is only one possible spin orientation in the lowest Landau level (1 = 0). 
With this result it is straightforward to calculate expectation values like 



2 In our convention three- vectors such as p = (p x ,Py,Pz) an d 7 = (7$, 7y,7^) are the 
contravariant components of the corresponding four-vector and thus have Lorentz indices 
% = 1, 2, 3 upstairs, i.e. p x = p 1 etc. We use the Minkowski metric diag(+, — , — , — ) so that 
Pi = —p l and 7i = —7* for i = 1,2, 3. 
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(*(a?)7 i T5*(«)> = ~ tr [iS{x,x)f^ & ] 



(A5) 



where the trace is over 7-matrices. Since 7^75 contains an odd number of 7-matrices, 
only the term in the integrand in Eq. ( |A4j ) , which is odd in p Q , can contribute. 
Evidently it is zero for a vanishing chemical potential, showing in a more formal 
way that the magnetization term of Semikoz and Valle || cannot be correct. 

It is often more convenient to label the Landau levels with an orbital quantum 
number n = 0, 1, 2 . . . and a spin quantum number A = ±1. The energies are then 
En,\,p z = m2 + vl + \qB\{2n + 1 — A). For a charged Dirac fermion / the net total 
number density (particles minus antiparticles) is 



N 



f-f 



\qB\ 
2tt 2 Jo 



poo 00 r 

J ^ n \ 1 1 



n=0 A=±l 

while the net number density in the lowest Landau level is 

\qB\ 



(A6) 
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iV /-/ 



2tt 2 Jo 



dPz /f ( E o,x,pz) - /f ( E 0,1, P z 



(A7) 



These results allow us to relate the local terms of the neutrino self-energy to the 
total charge density or to the charge density in the lowest Landau level, leading to 
Eqs. i) and (|). 

For the non-local neutrino self-energy term it is convenient to start from the 
electron propagator in the Schwinger proper-time form, which can be written as 

0,0: 



iS(x', x" 



fx', x") 



-iplx 1 — x") 



iS(p) 



(A8) 



where (f)(x',x") is a gauge-dependent phase factor. The gauge-independent and 
translationally invariant part of S is 



iS(p) = iS vl ,c(p) - / F (po) iS vac (p) - iS* ac (p) 



(A9) 



where 



roo e iqBscr z 
^vac(p) = / ds — — 

Jo cos{qBs) 



exp 



is p\ 



tan(qBs) 2 2 N 

p, — m + is 

qBs 



-iqBsaz 



1P\\ 



cos(qBs 



■ 7pj_ + m 



(A10) 



where for general four- vectors a and b, a ■ b\\ = a bo — {B ■ a)(B ■ b) and a ■ b± = 
a ■ b — (B ■ a)(B ■ b). The real combination that occurs in the thermal part of 
Eq. ( |^9| ) is obtained by extending the s-integral in Eq. ( A10 ) from —00 to +00. 
In the integrand of Eq. ( |A10|) there are poles and essential singularities on the real 
s-axis. They have to be avoided by taking the integration contour in the lower half- 
plane for positive s (see e.g. Ref. |L7j for a discussion of this contour). Therefore, to 
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get a real quantity for the thermal part, this contour has to go in the lower half-plane 
for negative s as well. 

The W boson propagator has a similar form but with a different tensor structure. 
In a closed loop, the gauge- dependent phase factors (f>(x, x') cancel and the result 
is explicitly translationally invariant. The contribution from thermal W bosons is 
Boltzmann, suppressed by a factor e~ mw ^ T and can be neglected. Expanding the 
W propagator in both the momentum transfer and the B field we find that the 
leading Independent 0(m^/)-term is local and that the first non-local 5-dependent 
term is 0{m^). The local term vanishes in a CP symmetric plasma. Therefore, 
when calculating the neutrino self-energy to order we may use the zero-field W 
propagator. 

The advantage with the Schwinger proper-time form over the Furry picture is 
that the gauge-dependent phase factor disappears automatically and we do not have 
to match the wave functions of the electron propagator (i.e. the Landau levels) with 
the ones of the W propagator in the zero field limit (i.e. plane waves). 

With the propagators in Eqs. (0) and ( |A9| ) it is possible to perform the loop 
integral over the three-momenta in Eq. (|2|) explicitly, but the result is still fairly 
complicated. It simplifies considerably in the linear-field approximation (B <C T 2 , 
B <C m 2 ), where we have, from the W^-e-loop: 



]bubble = " 2^ / MPo) £L dS e " (p2_m2 V [TP + m + iqB ■ as( 1Pll + 



x 



9nv{k~p) — (k — p) M (& —p)i 



(All) 



After adding the Z-i/-loop and keeping only the leading high-temperature piece, 
we obtain the result in Eq. (12). However, Eq. ( All ) is valid also for temperatures 



lower than the electron mass m. It contains corrections to Eq. ([12]) , which can be 
important if T < m. 
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